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Introduction. One-dimensional quantum gas has been 
attracting considerable attention since its seminal exper- 
imental demonstrations by using cold atoms [HQ- To 
date, one of the bosonic gas models extensively studied 
is the Lieb-Liniger model 0], which is solvable and char- 
acterized by a tunable coupling constant. 

In the case of repulsive interaction, much literature is 
found both for the zero-teinperature region and the 
finite-temperature region |7l-[ig|. The Yang- Yang equa- 
tions 0, whose solutions describe the thermodynami- 
cal behavior of the system, have been demonstrated ex- 
perimentally Il0| . Generalization to anyonic systems is 
also found in |ll| - [l3| . In parallel, the case of attrac- 
tive interaction has been investigated [l3 - [2l| in the last 
decade, which uncovered various exotic features includ- 
ing bound states with complex quasimomentum [ij- 
IITI [2H . Phase structure with respect to the coupling 
strength is found in the paradigm of Gross-Pitaevskii 
mean- field theory [l^ [2^ . 

In this Brief Report, we present two simple but no- 
table results about the Lieb-Liniger model valid for weak 
attractive interactions. One of them is the explicit form 
of the quasimomentum and energy spectrum, and the 
other is the equation of state in the high-temperature re- 
gim i which resembles the van der Waals equation of state 

Energy Spectrum in the Weak Coupling Limit. The 
Lieb-Liniger model 0] describes an A^-partite bosonic 
system in one-dimensional space with point interactions. 



It is governed by the Hamiltonian 



H 



2m ^ dx^ 

z— 1 * 



c'^5{xi 



Xj), 
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where m is the mass of the particle, c > is the 
coupling constant of the attractive 6 interaction, and 
the variable Xi G [0, L] represents the coordinate of 
the ith particle. The eigenfunctions iJj are symmetric 
• ■ ,Xi,--- ■ ■) = ■ ■ ,Xj,--- ,Xi,---) and obey 
the periodicity conditions ■ ■ ,Xi, - ■ ■ ) ~ ■ ■ ,Xi-\- 
L, • ■ • ) for all i. For convenience we hereafter work with 
the unit h = m = 1. 

Following the standard treatment of the model, we 
adopt the Bcthe ansatz. 



(2) 



for the simplex region 0<2:i < X2 < •■• < xn < L. Here 
ki is a quasi-momcntum (rapidity) and the summation 
is over all permutations a of the particles i t-^ a{i) on 
which the coefficients depend. Plugging ^ into the 
Schrodinger equation Hip = Eip, we obtain the Bethe 
equations [1] 



Lk^ = (2nj + N - 



1 ) TT + 2 arctan Si, 



(3) 



Here, rii is an integer, Sij is the rescaled relative rapidity 
defined by 
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{ki kj^ / C 
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for c ^ 0, and the arctangent takes the principal value, 
jarctani^i jl < 7r/2. The energy E then reads E = 



2 



In what follows, we consider the integer set {ui} so that 
ki G R for all i for the (dimcnsionless) weak coupling, 

e := cL/tt < 1. (5) 

We can confine ourselves to the case of the ordering, 

fci < fc2 < • • ■ < fcAT, (6) 

by relabeling the indices i appropriately. Note that (|6]) 
implies 6i,j > for i > j. We also assume that in the 
weak regime kt are all regular with respect to c. 

Two remarks are in order. First, there actually exist 
{rii} for few-partite systems on which the condition ([5]) 
is satisfied. For example, it is shown in Q that fc^ 6 M 
for bipartite systems for e < A/n. For tripartite systems, 
the same holds in the weak-coupling regime if the system 
has nonzero {rii} at the noninteracting limit c — ?■ fl5l |. 
Second, under the regularity assumption the ordering © 
implies no level crossing in the weak regime. 

We now show the following proposition. 

Proposition 1 For e <^ 1, the reseated relative rapidity 
and the coupling have the tradeoff relation, 



£(5, J = 2(n, ~ fij) + 0{e) 



(7) 



with 



Hi :~ Hi + i — 1. (8) 
Proof. To show this, we first note that for c — ?► 0, 

J tt/2 for i > j. 



arctan S 



1,3 



— 7r/2 for i < j. 
Substituting this to the Bethe equations ([3]), we find 



(9) 



, 27r_ 



(10) 



with rii given by ([5]) . This implies that the number rii is 
nothing but the quantum number of the free case c = 0. 
Since ^ holds for c = 0, for the weak coupling regime 
the regularity condition assures that 



and also that 



TT-l < 712 < ■ • • < riN, 



27r 1 
h = —rii + ;^C'(e), 



(11) 



(12) 



on account of 0{c) = {l/L)0{e) from ©. We then ob- 
tain 



e5.,+i., = -{h+i - ki) = 2(71, - 77,) + 0(e), (13) 

TT 



which completes the proof. 



The tradeoff relation ([7]) implies \5ij\ = 2(77; — nj)/e + 
0(1) ^ 1 or 1/ \5i.j \ = C(e), which is useful to approxi- 
mate ki. Indeed, with 



j-l/.T + 71/2 + 0(17x3), .T>1, 

arctana; = < ' ' ) ' ' (14) 

l-l/a--7r/2-t-0(l/a;3), x < -1, ^ ' 



substituted for Eq. ([3]), we find 



2^ 2^1 1 



Utilizing ([7]) again, we can eliminate 5ij from (|15p to 
arrive at the following theorem. 

Theorem 1 For e 1, the rapidity is approximated by 

This is one of the announced results of this paper. 

Let us now derive the approximated energy spectrum 
E to the order of 0(e). Summing up all the squared 
rapidities (fT6|) . we obtain 



ni,n2,....nN 



^ n 1-2 



z— 1 I— 1 'j^i 

^ 2^2^2 N{N-1) 



(17) 
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where we have used 

N _ N _ N 

rii \ - \ - 7li \ - \ - 77, 



77; — 77, -^-^ -^-^ 77,; — 77, -^-^ 77, — 77^ 

w _ _ 
^ — ^ ^ — ^ — 

= N{N-\). (18) 

The energy spectrum (fTT]) is additive for the particles, 
and we may rewrite it as En^^n^_...nj 



YliLiEn, with 



Ef: 



27r2 _2 - 1 



(19) 



The system thus behaves as an assembly of non- 
interacting particles of the energy En, consisting of the 
kinetic part 27r^77^/L^ and the averaged potential part 
—c{N — 1) / L which is proportional to the number den- 
sity {N - l)/L. 
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FIG. 1. (Color online) Pressure of one-dimensional Cs gas 
as a function of the one-dimensional scattering length am, 
evaluated for N = 25, L ^ 30 /im, T = 0.16 uK. For lower 
values of am, we find that P may become negative, which 
occurs before reaching the shaded range am < 2 x 10 fim 
where the weak coupling condition ((SJ is no longer valid. 

Equation of State. Let us consider the equation of 
state of the system for the high-temperature region: 

P :== l/ksT < L^/2tt'\ (20) 

where fee is the Boltzmann constant and T is the tem- 
perature. Under the Maxwell-Boltzmann distribution, 
the partition function Z becomes 

N 

tN a\-N/2 PcN'^ 

(27r/5) ' exp . (21) 

L 

To obtain this, we approximated the sum over {n;} by the 
integral over pi := iirni/ L, ignoring the case of complex 
rapidities (i.e., pi = pj for some which is measure 
zero. 

Since the volume of the system is given by L in our 
one dimensional model, the pressure P is evaluated as 
P = ^-jj-lnZ. We then obtain the following theorem. 

Theorem 2 The equation of state in the high- 
temperature region reads 

(P + CjjjL^NkBT. (22) 

This is in fact the van der Waals equation without the 
volume correction. The absence of the volume correction 
is understood by the fact that the Lieb-Liniger model 
adopts only point particles. 

The one-dimensional system of bosons with attractive 
interactions of our interest has actually been realized as 
the super- Tonks-Girardeau (sTG) phase of an ultracold 
gas of cesium (Cs) atoms using the confinement-induced 
resonance In this experiment, one observes a multi- 
ple of elongated tubes, each consisting of approximately 
25 atoms along the size of 3 x 10 ^m on average. To an- 
alyze this in our context, we take these numbers for N 



and L, respectively, with the mass m « 2.2 x 10 kg. 
From (|20p . we then find that our equation of state ([^^ 
is valid for T > 2n^h^/mL^kB « 8 x 10 pK. On the 
other hand, according to [IJl our coupling c is related to 
the tunable one-dimensional scattering length oid > 
by c = 2h^ /maiY). Thus our weak coupling condition ([5]) 
implies aiu ^ 2L/7r « 2 x 10 /^m. In contrast, we find 
fliD ~ 3 X 10^ nm in the experiment j23| . which is way 
out of the valid range. 

If we are allowed to extend the range by tuning the 
parameters of experiment such as the strength of applied 
magnetic field appropriately without spoiling the sTG 
phase up to the valid range, then our formula suggests 
an intriguing possibility. This is seen in the behavior 
of pressure P (see Fig. [T]) which becomes negative for 
lower values of aiD. This takes place when the pressure 
correction term surmounts the ideal gas term in (|22p . 
suggesting a phase transition there. This may also be 
realized, e.g., by increasing the number density N/L even 
if aiD is fixed. 

Conclusion and Discussions. In this paper we ob- 
tained the approximated energy spectrum of the one- 
dimensional attractive Bose gas where particles interact 
weakly with each other via the contact S potential. The 
coupling e and the rescaled relative rapidity Sij fulfill 
the tradeoff relation (Proposition [T]). As a result, we ob- 
tained the closed expression of the approximated rapidity 
(Theorem [T|). Based on this, we found that the equation 
of state is the van der Waals equation without the volume 
correction (Theorem [2]). 

The induced correction term is proportional to the 
attractive coupling c with the squared number density, 
which is precisely the same as the standard correction. 
This correction is conventionally ascribed to the average 
interparticle potentials in the classical thermodynamics 
derivation, whereas in our quantum-mechanical model 
the effect is taken care of by the nontrivial boundary con- 
ditions of the S potential imposed at the positions of the 
particles. The valid range of our analysis does not seem 
to overlap with the parameter range of the present exper- 
iment with ultracold Cs gas, but if it can be extended we 
may observe a novel phase transition as a consequence of 
the quantum boundary effect. 

In closing, we mention that the 5 potential is not the 
only possible point interactions admitted in quantum me- 
chanics. Indeed, it has been known that in one dimension 
we have a U(2) family of point interactions each charac- 
terized by distinct boundary conditions [25j including the 
S potential as a special case. Despite that all of them are 
equally zero-range, these potentials can in principle give 
rise to different spectra through the non-trivial boundary 
conditions, bearing a variety of palpable physical effects 
as pointed out in j26l - l29l | . We may therefore expect sim- 
ilar exotic outcomes to appear, for instance, in the equa- 
tion of state when such potentials other than the S are 
considered. 
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